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ABSTRACT 

We report on a detailed study of the water maser kinematics and radio contin- 
uum emission toward the most massive and young object in the star-forming region 
AFGL 2591. Our analysis shows at least two spatial scales of multiple star formation, 
one projected across 0.1 pc on the sky and another one at about 2000 AU from a ZAMS 
star of about 38 Mq. This young stellar object drives a powerful jet- and wind-driven 
outflow system with the water masers associated to the outflow walls, previously de- 
tected as a limb-brightened cavity in the NIR band. At about 1300 AU to the north of 
this object a younger protostar drives two bow shocks, outlined by arc-like water maser 
emission, at 200 AU either side of the source. We have traced the velocity profile of 
the gas that expands along these arc-like maser structures and compared it with the 
jet-driven outflow model. This analysis suggests that the ambient medium around the 
northern protostar is swept up by a jet-driven shock (> 66 km s~^) and perhaps a 
lower-velocity (~ 10 km s~^) wind with an opening angle of about 20° from the jet axis. 

Subject headings: stars: formation — stars: individual: AFGL 2591 — ISM: kinematics 
and dynamics — masers — techniques: high angular resolution 



1. Introduction 



The study of high-mass ( > 8 Mp)) star forrn ation still faces fundamental questions such as 



(e.g., IZinnecker &: Yorke 



2007; 



Beuther et al.ll2007l ) what is the origin of the mass reservoir for the 



protostar (whether the protostar can gather mass only from its own cocoon or competes with the 
nearby protostellar companions to accrete mass from the whole molecular clump) and what is the 
accretion mechanism (whether the mass is conveyed to the protostar through a stable accreting 
disk or via episodic accreting flows while the protostar crosses the denser portions of the molecular 
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clump). To examine these questions, observations on linear scales matching the sizes of the molec- 
ular clumps (0.1-1 pc) need to be combined with a closer view of the protostellar surroundings 
(the hot molecular cores, HMCs, with sizes of a few thousands of AU) to establish the presence of 
multiple sources, their kinematical structures, and a possible overall interaction. 

Focusing on a single protostellar core. Very Long Baseline Interferometry (VLBI) proper motion 
studies of maser spots in conjunction with sub-arcsecond imaging of the free- free emission of very 
young H li regions and thermal jets have proven a powerful tool to investigate the kinematics and 
the structure of the gas close to the protostar(s). For instance, our VLBI multi-epoch observations 
toward the massive young stellar object (YSO) G23.01-0.41 have shown three different kinematical 
structures traced by the H2O, OH, and CH3OH masing molecules within 2000 AU from the central 
YSO. While the water masers trace a bipolar flow associated with a thermal jet, the hydroxyl 
masers belong to an expanding layer in front of the central source. Furthermore, the methanol 
masers appear to undergo both rotati on and expa.nsion i nside a toroidal structure traced in NH3 
lines and extended over about 0.1 pc (jSanna et al.ll2010bl ). 



The present paper focuses on the water masers and radio continuum toward the high-mass star- 
forming region AFGL 2591. The distance to this region has been substantially underestimated in 
the past and we have recently measured an accurate value of 3.3 kpc with the trigonoi netric parallax 
of its 22.2 GHz H9.O inasers (Rygl et al. 2011, submitted). Based on the data of iMinh fc Yang 
(j2008l ) and iLada et al.l (119841 ). we find a revised value for the large-scale (< 10 pc) clump mass 
of about 2 X 10^ and a bolometric luminosity inferred from IR data of 2 x 10^ Lq. From 
observations of the ^^CO rovibrational lines, the sys temic velocity of the reg ion (Vgys) with respect 
to the local standard of rest (LSR) is —5.7 km s~^ ( van der Tak et al.lll999l ). 



AFGL 2591 sho ws energetic star formation activity with hot core ernission detected in several 



mole cular lines (e.g.. Ivan der Tak et al.lll999l : lDoty et al 



2002 



Benz et al 



2011 



20071 : 



van der Wie^ 



. A number of compact radio continuum sour ces were resolved in the region (jCampbel 



Trinidad et al 



masers (jTofani et al.lll995l : 



Trinidad et al.l 



2OO3I ) . which establishes the occurrence of multiple star 



et al ■ 
19841 : 



2003 : Ivan der Tak &: Menten 20051 ) and associated with distinct clusters of water 



formation. AFGL 2591 exhibits a powerful CO outfiow, which is extended more t han 5^ x 5^ 



(^ 5pcx5pc), with the blue lobe toward the west and the red lobe to the northeast (e.g.. lLada et al 



1984 



Hasegawa fc Mitchelllll995r). This qutflq w activity is also associ ated with several Herbig-Har o 



objects toward the west (jPoetzel et al.l Il992l ) and H2 bow shocks (jTamura &: Yamashital Il992l ). 
marking the interaction between outflowing material and the surrounding molecular envelope. 
There is evidence for an outflow cavity that enhances the chemistry along the outflow walls far 
from the dust peak (e.g.. iBruderer et al.ll2009l ). This outflow cavity is obser ved as NIR loops , asso- 
ciated with the blueshifted outflow lobe, at the h igh resolution of 170 mas by lPreibisch et al.l (j2003l . 
cf. also iTamura et al.lE99ll : iMinchin et al.lll99ll ). The apex of these NI R loops is centered o n the 
most compact radio continuum source in the region called VLA-3 (after iTrinidad et al.ll2003l ). 



The observations reported here focus on the water maser emission within 2" from VLA-3. In 
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Sect. 2, we describe our Very Long Baseline Array (VLBA) observations of the 22.2 GHz H2O 
masers together with the archival Very Large Array (VLA) A-configuration observations of the 
radio continuum emission at 1.3 and 0.7 cm. In Sect. 3, we illustrate the spatial morphology and 
kinematics of the maser emission and present results from archival VLA observations, constraining 
the properties of the radio continuum associated with the masers. Section 4 discusses the properties 
of each cluster of water maser emission detected in the region and draws an overall view for the 
purposes of star formation. The main conclusions are summarized in Sect. 5. 



2. Observations and Data Analysis 

2.1. VLBA observations: 22.2 GHz H2O masers 

We conducted VLBaI observations in the K band to study the 616 — 523 H2O maser emission 
(rest frequency 22.235079 GHz) toward AFGL 2591. The VLBA observations were scheduled under 
program BM272H at four epochs: 2008 November 10, 2009 May 6, 13, and November 13. While 
these observations were optimized to measure the annual parallax and Galactic proper motion of 
the source (Rygl et al. 2011, submitted), in this paper we focus on the internal kinematics of the 
water maser emission. 

We performed phase-referencing observations by fast switching between the maser target and 
the ICRF calibrator J2007+4029. That allowed us to determine the water masers absolute posi- 
tion within an accuracy of ±1 mas. Two fringe finders (3C345 and 3C454.3) were observed for 
bandpass, single-band delay, and instrumental phase-offset calibration. Further details about these 
observations as well as the distance measurement to AFGL 2591 can be found in Rygl et al. (2011, 
submitted). We employed four adjacent IFs of 8 MHz bandwidth in dual circular polarization, each 
one split into 256 spectral channels. This receiver setup provided both a large-enough bandwidth to 
increase the SNR of the weak continuum calibrators and a 0.42 km s~^ channel width to adequately 
sample the maser lines. The third IF centered on an LSR veloci ty (Vt,sr.) of -5 - kni s~^ covered 



the previous detections of water maser emission in the region (jTrinidad et al.M2003l ). The data 
were processed with the VLBA correlation facility in Socorro (New Mexico) using an averaging 
time of about 0.9 s which limited the instantaneous field of view of the interferometer to about 2" 
(i.e. without significant amplitude losses). Data were reduced w ith the NRAO Astr onomical Image 



Processing System (AIPS) following the procedure described in lReid et al.l (I2009a|). 



The natural CLEAN beam was an elliptical Gaussian with an FWHM size of 0.80 mas x 
0.40 mas at a P. A. of -14.9° (east of north), with small variations from epoch to epoch. In each 
observing epoch, with an on-source integration time of about 2 h, the effective rms noise level on 
the channel maps was about 0.01 Jy beam^^. The total-power spectrum of the 22.2 GHz masers 



^The VLBA is operated by the National Radio Astronomy Observatory (NRAO). The NRAO is a facihty of the 
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toward AFGL 2591 is shown in Figured! 



2.2. VLA archival data 



We retrieved VLA A-configuration observations toward AFGL 2591 at 1.3 and 0.7 cm from 
the VLaI^ archive for the purposes of a detailed comparison with the high-resolu tion VLB A maser 
dat a. These observations were s cheduled on 1999 June 29 (jTrinidad et al.ll2003l ) and 2002 March 
23 (jvan der Tak &: Mentenll2005l ). respectively. We recalibrated these data with the NRAO AlPS 
software package using standard proce dures and cross-calibr ation between the strong H2O maser 
lines and the continuum at 1.3 cm (e.g.. lReid &: Mentenlll990l ). For the 1.3 cm data, we cleaned the 



VLA maps with a circular restoring beam, equal to the natural-beam minor axis for a "ROBUST 
0" weighting (Tabled]). 



3. Results 

3.1. 22.2 GHz H2O masers 

We imaged a range of LSR velocities between -40 to -|-25 km s^^ with a field-of-view of 
l'.'4 X 1^.^4 about t h e VL A-3 radio source. These limits include the emission previously reported 



by iTrinidad et al.l (120031 . their Table 3 and 4) toward VLA-3. In the following, the term feature 



refers to spots spa tially overlapping i n contiguous velocity channels, that we treat as an individual 



masing cloud (see lSanna et al.ll2010al for a detailed discussion). 



We identified 80 distinct water maser features distributed within an area of about 0'.'77 x 0'.'48 
(~ 2500 AU X 1500 AU). They are mainly grouped in five clusters labeled in Table [2] north-west 
(NW), middle- west (MW), north-east (NE), middle-east (ME), and south-east (SE) in Figure [2li. 
The properties of individual features are presented in Table [2] and show peak brightness ranging 
from 22 Jy beam~^ to a detection limit of 0.05 Jy beam~^ (5o"). Water maser features redshifted 
with respect to the Ygys (of —5.7 km s""*^) belong to the MW cluster, whereas the blueshifted 
emission is clustered toward the SE. The full spread in LSR velocities ranges from —0.4 km s~^ 
for the most redshifted feature (# 13) to —34 km s~^ for the most blueshifted one (# 77). The 
angular distribution and full-space motions of the water masers (i.e. line-of-sight velocities plus 
proper motions) are plotted in Figure [2^. With a time baseline of 1 yr, we measured relative proper 
motions of individual features with an uncertainty less than 3 km s^^. Positions and velocities in 
Table [2] are relative to the compact feature ^ 4 which was stable in the spectral domain over the 1 yr 
observations. The magnitude of relative proper motions ranges from 5.0 km s~^ for feature 7^ 17 
to 56 km s~^ for feature # 54. Relying on the symmetry of the maser spatial distribution, we 
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have corrected the proper motions by the average velocity of all features of Vx = 8.8 it 0.4 and 
Vy = 9.0 ± 0.4, to approximate the actual motions with respect to the center of expansion. 

In Figure [2)3, [2}3, and [2]f, we zoom in on the three water maser clusters labeled MW, NE, 
and SE, showing a symmetric expansion away from their centroids. The NE and SE clusters are 
particularly interesting due to their full kinematics and are discussed in detail in Section [4 . 1 1 and Hr2| 
respectively. 

In Figure [2jl-e, we overlay the radio continuum emission described in Section 13.21 with the 
water maser emission from the second epoch observations. This maser emission was summed over 
the full velocity range of the SE cluster (from Vlsr of —34 to —14 km s~^). Its projected shape 
along the plane of the sky delineates a straight "V" (hereafter the II2O V-shape) with a full opening 
angle of 110° it 10° and a P.A. of 256°. The variation of LSR velocities along this V-shape is quite 
regular, ranging from a maximum of —13 km s~^ (feature ^ 78) at the northern side of the V-shape 
to a minimum of about —33 km (feature # 77) toward the southern side. 



3.2. Radio continuum emission 



Our measurements of the 1.3 and 0.7 cm continuum emissions are summarized in Table [T] 
and plotted in Figure [2]i-e. Published VLA C-configuration observations at 0.7 cm show that 
significant emission is resolved out with the VLA A-configuration, missin g more than half of the 
flux measured in the compact configuration (Ivan der Tak &: Mentenll2005l ). The relative position 
of the 1.3 cm continuum with respect to the reference maser feature in Table [2] is obtained by 
matching the overall water maser distribution between the VLB A and VLA-A maps. We note that 
the NS spread of the maser emission overl aid on the radio cont inuum (~ 50 mas) have remained 
almost stationary over about ten years (cf. iTrinidad et al.l 120031 ). 



The peak position of the 0.7 cm continuum is assumed to be aligned with the 1.3 cm peak. 
The three-lobed geometry of the 7 mm emission follows the 1.3 cm elongation toward the east 
and it branches off toward the west matching the H2O V-shape (Table [I]). While this matching is 
particularly striking, we note that more sensitive, Q-band, observations are mandatory to confirm 
the spatial structure of the 0.7 cm emission. Because of the stability of the water maser distribution, 
but also the similarities between the morphology of the water maser emission and the 0.7 cm 
continuum (Figure [2^), we estimate a relative position accuracy between the radio continuum and 
water masers of itO'.'Ol. 

Previous observations spanning a broad-band spectrum between 6 cm to 1.3 mm show that 
the spectral index (S oc u"') for wavelengths longer than 7 mm is dominated by free-free emission 
(7 ^ 0.9), whereas at shorte r wavelengths the flux density is mostly due to dust emission (7 > 2; 
van der Tak &: Mentenll2005l and reference therein). 



- 6 - 



Discussion 



4.1. Outflow emission in the NE H2O cluster 

In this section, we discuss two basic characteristics observed in the NE water maser cluster 
(Figure [2)3) . Firstly, we detected two arc-like filaments of water maser emission elongated in the 
NE-SW direction with no emission in between them. Secondly, these arc-like filaments expand from 
each other with proper motions increasingly aligned with the NE-SW direction as one approaches 
this direction. Bipolar distributions of water masers wit h expanding prope r motions have been 



reported in a large nu mber of star- forming regions (e.g., iGwinn et al 



1992 



Sanna et al. 2010b 



Moscadelli et al.l 120111 ). These detections were interpreted as footprints of the interaction between 
collimated jets and/or wide-angle winds from YSOs with the surrounding molecular envelope. 
Following this evidence, we want to show that the overall kinematics of the NE water maser cluster 
can be explained by assuming that a single object, in between the maser emission, emits a high- 
velocity, bipolar jet along the NE-SW direction (Figure [2)3 and [3^). At first, we will compare 



( 


2001 


) and 


Lee et al. 


( 


2001) 



the interaction between a bow shock and the ambient material gives an upper and lower limit to 
the velocity of the processed material. This processed material is a mixture of material that was 
already part of the bow shock and material incorporated from the ambient medium as the shock 
front expands. The shock front delineates an arc-like structure (as do the masers) and the velocity 
component along the jet axis d ecreases roughly as r ~^/'^, where r is the offset distance from the 



head of the jet (e.g.. Figure 2 in lOstriker et al.ll200ll ) 



The similarities of the spatial elongation and velocities of the water masers would be naturally 
explained if the NE-SW direction is associated with an outflow in this direction. In Figure [3^, we 
define an elliptical symmetry that matches to first order both the width and the NE-SW elongation 
of the water maser distribution with a simple central symmetry (major axis and eccentricity on top 
of Figure [3ti.) ■ This elliptical pattern allows us to define a new reference frame with origin at the 
center of the ellipse (Ax = 625 mas and Ay = 260 mas w.r.t. the reference feature in Table [2]) and 
x, y axes along the minor (R) and major (z) axes of the ellipse, respectively. This reference frame 
is tilted by 22° to the east of the north direction. Due to the small dispersion of LSR velocities 
(1.5 km s~^) centered at about —7.9 km s~^, we also assume that all the maser features lie along 
the same plane. We want to study the velocity profile of the shock fronts traced by the water maser 
proper motion s with respect to the ne w z-R reference f rame. In the following, we use the same 
terminology of lOstriker et al.l (j200ll ) and lLee et al.l (|200ll ) in order to facilitate the comparison with 
their jet model. 

In Figure [3t), we report the behavior of the velocity components (Vz) of the water masers 
along the major axis of the ellipse (z). The idea that the z-direction is a preferred outflow direction 
is supported by the velocity proflle in Figure [8)3, which resembles that expected for a ballistic 



bow-shock produced by a jet along the z-direction (cf. Figure 4 of iLee et al.ll200ll ). For such a 
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bow-shock model velocities are constrained in Figure [8)3 by two solid lines and the expected gas 
velocity should lie b etween them. These limits are derived from equations (22), (18), and (20) of 
Ostriker et al.l (|200ll ) and define the bow shock profile z = f{R), the mean velocity of the mixing 
of jet and ambient material v'^ (solid line) in a shell, and the velocity of newly shocked ambient 
material u'^ (dashed steeper line), respectively. These equations are parameterized by the jet radius 
Rj, the velocity of the bow shock Vg, and the ejection velocity /3cs of the material processed b y the 



working surface, expressed as a multiple of the isothermal sound speed (jOstriker et al 



2001 



For 

the radius of the high- velocity "core" of the jet (Rj), we set the minimum projected distance to the 
R-axis measured with the masers of 1.8 mas (~ 6 AU). This value implies a width-to-length ratio of 
the jet from the center of the ellipse of 0.05, consis tent with that found i n the literature for H2 jets 
at typical scales of 10^ - 10"^ AU (e.g., HH 212 in lzinnecker et ahlEoQa ). We stress that Rj is the 
only constraint we put on the model, whereas we consider the b est value of the ratio /3cs/vs = 0.5 
and /3 = 4.1 inferred from the simulations by iLee et al.l (l200lh. By setting th e tip of the bow 

we find a good 



Ostriker et al 



2001 



shock at the maximum maser projection along z (the "0" in 
agreement of the measured maser velocities (black squares in Figure [8)3) with the jet-driven outflow 
model. We note that the maser ve locities ar e close r to the ambient material velocity (dashed line) 
in agreement with simulations by iLee et al.l (120011 . their Figure 4). This behavior may reflect an 
incomplete mixing of material already in the bow shock shell and material added from the ambient 
medium. The /3cs/vs ratio implies a bo w shock speed of abou t 66 km s~^ for a sound speed of 8 km 
at the cooling cutoff of 10^ K (see Ostriker et al.|[200ll ). On the contrary, in the alternative 
scena rio of a wide-an gle wind the Vz motions would increase linearly with z, that is not observed 
(e.g., iLee et al.ll200ll their Figure 9). This suggests that the maser cloudlets are dragged along the 
z-direction by a jet-driven shock component. 

On the other hand, the jet-driven outflow model predicts an upper limit to the transverse 
velocity component Vr, with respect to the jet direction of the order of the sound speed (~ 8 km 
s~ , see Figure 3 in lOstriker et al.ll200ll and Figure 4 in iLee et al.l[200ll ). This limit is intimately 
related to the assumption that the transverse radial forces are dominated by thermal pressure 
gradients that act at the head of the jet. In Figure [3]3, we measure the greatest deviations from 
the jet model far from the head of the jet, at an inclination of about 20° from the jet axis. This 
behavior may suggest that, away from the head of the jet, thermal pressure gradients no longer 
dominate the expansion. One possibility would be that a low-velocity (< 10 km s~^) wind with a 
wide opening angle (ca. 20°) acts with the jet simultaneously. In the wind-driven outflow model, 
the expansi on of the outflow away from the central axis is driven by the ram pressure of a wide- 
angle wind (IShu et al.lll99ll ): the Vr velocity is ni aximum far from the tip of the outflow shell and 



goes to zero at its tip (Figure 9 in lLee et al.ll200ll ). The simultaneous presence of the two types of 
outflow mechanisms would agree with the hypothesis that any outflow contains a wide-angle wind 
at some level, in order to e xplain typical wi dth-to-length ratio of about 1 : 10, much higher than 
in a "pure jet model" (e.g.. lArce et al.ll2007l ). 



Finally, if the single jet-driven model is correct, we can also get an estimate of the dynamical 
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age of the jet-event exciting the water masers. The velocity of the jet (vj) is related to the bow shock 
speed (vg) by the ratio of the je t-to-ambient density , where Vg approximates vj if the jet density 
exceeds the ambient density (see Ostriker et al.|[200ll ). By considering the lower limit vj ~ Vg, the 



inferred value of Vg 
dynamical age of t = 

estiinates based on the expa nding motions of masers around forming YSOs (e.g., iTorrelles et al, 



= 66 km s , and the value of the higher maser projection along z, we obtain a 
Zmaj/vj — 14 yr. On the one hand, this dynamical age is consis tent with other 



2003, 



Moscadelli et al.l 120071 ). On the other hand, given how unlikely it would be to observe such 



a short-lived event, it strongly suggests that we are tracing a recurrent phenomenon of ejection of 
matter from an YSO. That would favor a pulsed-jet paradigm against a steady-jet one as supported, 
for ii istance, by the recurre nt knots and bow shocks observed in the prototypical Herbig-Haro object 
212 (jZinnecker et al.lll998l ). This interpretation naturally explains also the eastern water masers in 
Figure [3^ (features ^ 38, 39 in Table [2]) and the masers at about 150 mas to the south of the NE 
cluster (the ME cluster in Table [2]). These masers would belong to an older jet emission from the 
same object. 



4.2. Outflow emission in the SE H2O cluster associated with VLA 3 



In this section, we discuss three characteristics observed in the SE water maser cluster: I. the 
projected V-shape of the water maser distribution on the plane of the sky, "pointing" toward the 
peak of the radio continuum emission (Figure [2]l-e); II. the water maser proper motions that mainly 
expand along the V-shape (Figure [2^); III. the line-of-sight velocity gradient observed through the 
V-shape, with velocities increasing to the south (Figure [2f)- In the following, we relate these 
properties with the large-scale outflow and the nature of the radio continuum emission and suggest 
that the H2O masers trace the edges of the outflow cavity associated with a massive YSO (see 
Figure H^) . 

The H2O emission from the SE cluster is blueshifted with respect to the overall velocity of 
the cloud (—5.7 km s~^) and covers a range of line-of-sight velocities cons istent with the small 



scale , blueshifted wing of the CO(3-2) outflow (between -10 and -40 km s ^ : iHasegawa &: Mitchell 



I995I ). The V-shape observed in the water maser emissi on opens towa,rd the west, where a limb- 



brightened cavity wa s previously outlined as NIR loops ([Tamura et al.l Il99ll : iMinchin et al.l Il991 



Preibisch et al.l l2003l ). These loops were explained as wind-bubbles, which were created as the 
blueshifte d (western) side of t he outflow pushed against the environment and swept up ambient 
material ([Preibisch et al.ll2003l ). We note that the wide opening angle of the NIR loops (> 100°) 
is consistent with that of the V-shape (110°). Actually, if the loops outlined the outflow cavity, 
one might expect the water masers along the outflow walls, where shocks are more likely. The II2O 
V-shape opens toward a position angle (P.A.=256°) in good agreement with the direction of the 
western side of the outflow, which has been detected in a number of tracers: i) the NIR bright 
loops, with a 262° < P. A. < 256°, surrounding the outflow cavity with a symmetry 
of 259° (e.g., IPreibisch et al]l2003l . see their Figure 3); ii) the blueshifted CO lobe with a P. A. of 
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about 260° (e.g.. iHasegawa &: Mitchell! Il995l. see their Figure 4); iii) the western, aligned, Herbig- 
Haro obje cts with a 258° < P.A < 2 61° (jPoetzel et al.1 Il992l ) and the hot H2 knots in the same 
direction (jTamura &: Yamashital Il992l ): iv) the elongated 1.3 cm continuum with a P.A. of 268° 
and a spectral index consistent with an ionized wind (Figure [2]1). If water masers would belong to 
the blueshifted outflow lobe, and might possibly amplify the continuum background at 1.3 cm, it 
would naturally explain the non-detection of redshifted maser emission, that would belong to the 
receding outflow lobe. Also, since the outflow axis is close to the line- of-sight (within 30° — 45°; 
Minchin et al.lll99ll : IHasegawa &: Mitchelllll995l : Ivan der Tak et al.lll999l ). this scenario agrees with 
the direction of the maser proper motions, that are mainly expanding parallel to the V-shape with 
a strong component toward the observer (Figure U]). The velocity gradient observed across the 
maser V-shape (Figure [2f) might be reproduced with a slow precession of the ou tflow axis about 
the mean direction Zout in Figure H] (as speculated at first bv lTrinidad et al.ll2003l ). 



The strong water masers and t he radio continuum emission from VL A-3 (Figure [2]i-e) suggest 
an association with a massive YSO (Ivan der Tak &: Mentenll2005l ) . If the free- free emission longward 
of 7 mm wavelength originates from an ionized wind, we can estimate a lower limit to the mass- loss 
rate (M^) of t he wind and put a constr aint to the the nature of the driving source. We apply 
equation (4) of lRodriguez Sz Cantd (jl983l ) and use the flux at 1.3 cm to estimate the mass-loss rate 
in the assumption of a spherical, isothermal, completely photoionized, uniform wind: 



10-5 M0 yr- 



1.7 



1.6 mjy 

3/2 



3/4 



22.2 GHz 



d 



3.3 kpc 



10^ km s- 



(1) 



where is the velocity of the wind. We have considered the accurate value of the recently 
measured distan ce to AFGL 2591 an d a wind velocity of several 10^ km s~^, reported i n bot h 
the [SII] lines bv IPoetzel et al.l (Il992l l and the IR ^^CO broad wings bv lvan der Tak et al.l \imd ). 
The assumption of a spherical wind is justified by the w ide opening an g le obs erved both with the 
maser emission and the NIR loops (> 100°). Note that, iTrinidad et al.l (j2003l ) used the lower flux 
at 3.6 cm to estimate the mass-loss rate, based on their model with a dusty disk that would have 
affected the flux at 1.3 cm. Still, at the new measured distance of 3.3 kpc, their modeled flux would 
decrease by a factor 10 and no longer be applicab le, while the spectral index fo r wavelengths longer 
than 0.7 cm suggests an origin in a ionized wind (jvan der Tak &: Mentenll2005l ) . The mass- loss rate 
from equation ([T]) implies a photoionizing ZAMS star of type 09, with a rate of ionizing photons 
(N^) of about 1.2 x 10^^ s~^ (jRodriguez &: Cantdll983l ). We compare this value with the estimates 
of mass and Lyman contii iuum flux derived by assuming a single star dominating the IR luminosity 
of the region. Fol l owing Ivan der Tak &: MentenI (|2005l ) , we used the theoretical HR diagram by 
Maeder &: Mevnetl (|l989l ) to estimate the mass and effective temperature of a star with an IR 
luminosity of 2 x 10^ L0; then, we used the stellar atmosphere models by ISchaerer &: de Koter 
(j 19971 ) to convert effective temperatures to Lyman fluxes. These new values of M and N/, are about 
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ss Mq and 2 x 10^^ s~^, respectively. The estimates of Lyman flux obtained independently from the 
radio continuum emission and the IR luminosity of the region agree within an order of magnitude. 
Since the former one is a lower limit derived close to the source whereas the later one gives an 
upper limit to the number of Lyman photons emitted over a larger field of view, this calculation 
provides evidence that a single mass ive object dominates the IR energetics of the region (according 
to the high-resolution NIR image of lPreibisch et al.ll2003l ). 



4.3. Clustered star formation around VLA 3 



We finally collect the information from the water masers and the radio continuum emission 
about the ongoing star formation in the region. Taking into account a possible bias of detecting 
only a small number of (bright) sources, there is evidence of at least two clustering scales for star 
formation in the region. One on a linear scale of about 0.1 pc on the sky, fr om continuum source s 
VLA-1, VLA-2, and VLA-3 detected at first bv lCampbell (|l984l . named bv lTrinidad et aLlbood ). 
The other scale is an order of magnitude smaller and comes from our observation of at least three 
centers of expansion that pinpoint two other centers of star formation (the NE and MW cluster), 
separated by about 2000 AU on the sky from the VLA-3 radio source. Note that, the NW maser 
cluster might be either associated with the MW cluster or trace a further young object. Our 
analysis in Section 14.21 shows that VLA-3 is the most massive (< 38 M0) and young object in 
the field dominating the IR luminosity. Several lines of arguments from the spatial scale of the 
maser distribution also strengthen the association between VLA-3 and the large-scale dynamics 
of the molecular envelope, as traced by both jet "footprints" (HH objects and H2 knots) and the 
wide-angle outflow seen at the sub-parsec scale (~ 0.5 pc from VLA-3). The other two objects 
indicated by water maser emission in Figure [2] (the NE and MW cluster) were not detected at cm 
and mm wavelengths. We give an upper limit to their flux density at 1.3 cm of about 0.5 mJy 
{3a). On the one hand, by assuming this radio continuum would come from a photoionized stellar 
wi nd, we can put an upper lim it on the spectral type of the exciting YSOs. Following equation (1) 
of [Rodriguez fc Cantol (119831 ) which gives the rate of ionizing photons required to fully ionize a 
constant velocity wind, we get a value of Nl ^ 10^^ s~^ that would imply a ZAMS star with 
a spectral type later than BO. On the other hand, the high, isotropic, maser-luminosity of the 
NE and MW objects (2.5 x 10~^ Lq and 3.3 x 10~^ L©, respectively) would suggest they are not 
even low-mass YSOs but possibly late B-type stars, by comparison with the SE cluster luminosity 
(1.2 X 10~^ Lp) an d the H2O luminosity usually associated wit h low-mass star s (< lO "'' Lq; e.g., 
Furuva et al.l l2003l and their erratumjj. In the formulation of iReid &: MoranI ( 19881 ). the maser 
luminosity is a consequence of the dissipation in the ambient cloud of the mechanical energy of the 
outflow from a young star, and thus could be related to the mass of the driving source. 



■^Note that from Trinidad et al. ( 20031 ) . the isotropic maser-luminosity associated with the VLA-2 radio source is 
of about 4 X 10~* L© (their Table 2) whereas VLA-1 does not show water maser emission. 



- 11 - 



5. Summary and Conclusions 

We have presented multi-epoch, VLBA, H2O maser observations toward the most compact 
radio source (i.e. VLA-3) in the high-mass star-forming region AFGL 2591. We have also compared 
the maser distribution with the brightness structure of the radio continuum at 1.3 and 0.7 cm 
observed with the VLA. Our main conclusions can be summarized as follows: 

1. We have detected three main clusters of water maser emission above a detection limit of 
0.05 Jy beam~^ {5a) within ~ 2000 AU on the sky from the VLA-3 radio source (the SE 
water maser cluster). Each water maser cluster shows internal expansion with velocities of a 
few tens of km s~^ (Figure [2^). A global view of the star formation activity in AFGL 2591 
shows at least two spatial scales of star formation, one projected across 0.1 pc on the sky and 
another one at a ten times smaller scale as indicated by water maser emission. 

2. The kinematics of the NE water maser cluster shows characteristics of two symmetric, jet- 
driven, bow shocks expanding about 200 AU from a still undetected protostar, with a velocity 
of about 66 km s~^ (Figure Eb). Our analysis of the water maser proper motions, compared 
to different outflow models, suggests that the primary wind may be a combination of a jet 
and a lower-velocity wind tilted about 20° from the jet axis. The short dynamical age (14 yr) 
of these expanding motions suggests they are probably related to recurrent ejection events 
from the central protostar (s). 

3. Our analysis of the radio continuum emission compared with the IR luminosity of the region 
shows that VLA-3 is likely a ZAMS star with a spectral type between 09-06 and a mass 
in the range 20 — 38 Mq. This young stellar object dominates the IR luminosity of the 
region. The water maser emission associated with VLA-3 expands toward the observer with 
velocities up to 50 km s~^ and traces a V-shape along the plane of the sky (Figure [2]i-f) . The 
orientation of the maser V-shape is aligned with the high- velocity wings of the CO outflow 
and the Herbig-Haro objects at a few tenth of pc from the star. The wide opening angle of 
the maser V-shape (~ 110°) agrees with that of the outflow as inferred from NIR emission 
along the edges of the outflow cavity. This water maser emission is likely amplified along the 
walls of the blueshifted outflow lobe from VLA-3 (Figure H]) . 



We thank J. M. Torrefies and A. P. Lobanov for helpful discussions in preparation. This work 
was partially funded by the ERG Advanced Investigator Grant GLOSTAR (247078). KLJR is 
funded by an ASI fellowship under contract number 1/005/07/1. 
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Fig. 1. — Total-power (Stokes I) spectrum of the 22 GHz H2O masers toward AFGL 2591 from 
the second (middle) epoch on 2009 May 6. This profile was produced by averaging the total-power 
spectra of all VLBA antennas, after weighting each spectrum with the antenna system temperature. 
The dotted line crossing the spectrum (at —5.7 km s~^) represents the systemic velocity (Vsys) of 
the molecular cloud hosting the star-forming region. The LSR velocity range from the SE water 
maser cluster is shown. 



- 15 - 



(b) - 



V 



7 0.65 e 

EAST OFFSET (arcaec) 




1300AU 





EAST OFFSET (; 



(a) 



■<1 



I 



20 29 24:90 



24?a6 24 84 

RIGHT ASCENSION (jaOOO) 




SE cluster 

- (f) ' 



0.63 
EAST OFFSET (. 



*LSH _^ 

km s 

-30,8 
-2B 
-25.2 
-22.4 
-19.6 
-16.8 
-14.1 
-11.3 
-8.5 
-5.1 
-4.5 
-3.9 
-3.3 
-2.7 
-2.1 
-1.5 
-0.9 
-0.3 



U V 



EAST OFFSET (ai 



Fie. Zi. 22.2 GHz H2O mascr kinematics and radio continuum toward the VLA— 3 source in AFGL 2591. Central Panel: Absolute positions 
(triangles) and internal proper motions of the 22 GHz H2O maser features (soc Tabic [2) . Colored cones are used to show both the direction 
and the uncertainty (cone aperture) of the 3— D velocities of maser features. The proper-motion amplitude scale and linear-size scale are given 
(10 mas — 33 AU). Different colors mark the maser LSR velocities according to the righthand side scale, centered on the systematic velocity of 
the HMC (green equal to -5.7 km s""'"). Upper Panels: Details of the proper motions in the NE and MW water maser clusters. Relative positions 
refer to the phase- reference feature # 4 in Table [21 Bottom Panels: from left to right, detail of the extended water maser emission (blue contours) 
overlapped to the VLA— 3 radio component of lTrinidad et al in the same reference frame of the upper panels. Blue contours are spaced 

by 3cr starting from a 3cr rms of 0.39 Jy beam""'". For the VLA— A radio continuum maps (black contours), plotted levels at 1.3 cm are spaced by 
Icr starting from a 3tT rms, whereas for the 7 mm map contours are -3cr, 3cr, 5tT, 7(7, and 9cr (see Table^J. A blow-up of the water maser proper 
motions associated with VLA— 3 is shown. 
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Fig. 3. — Symmetry and velocity profiles of the NE water maser cluster. Left Panel: Spatial 
distribution of water maser features (triangles) superposed to the elliptical toy model used to 
define the new z-R reference frame. The major axis and eccentricity of the ellipse are reported 
on top of the panel. Middle Panel: Velocity profile of the water maser proper motions along the 
z-direction in astronomical units. T he lines represe i it the upper (solid line) and lower (dashed 
line) limits to the bow shock model of lOstriker et al.l (j200ll ) discussed in Section [4.11 Right Panel: 
Velocity profile of the transverse component, with respect the the z-direction, of the water maser 
proper motions as a function of z. The lines have the same meaning as in the Middle Panel. 
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Fig. 4. — Sketch of the blueshifted jet/outflow system in AFGL 2591 emanating from the radio 
source VLA-3 and associated to the SE water maser cluster (Figure [2]l-f ) . The inclination of the 
outflow axis about the line-of-sight (~ 30 — 45°) and on the plane of the sky (~ 10°) is reported. 
The NIR loops are shown as a cross section of the outflow projected along the plane of the sky. The 
water maser emission associated with VLA-3 emerges from the outflow walls and mainly expands 
toward the observer (Section 14. 2p . 
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Table 1. Radio continuum toward the VLA-3 component in AFLG 2591 



Peak position 

Telescope A HPBW Image rms Ax Ay Fpeak Deconvolution 

(cm) mas X mas at ° (mjy beam~^) (mas) (mas) (mJy beam"-"^) mas X mas at ° 



VLA-A 1.3 80 X 80 0.19 633 ± 10 -145 ± 10 1.24 88 x 38 at 88° 

VLA-A 0.7 45 x 38 at - 34° 0.11 633 ± 10 -145 ± 10 0.95 65 mas for slice at 75° 

48 mas for slice at —48° 
62 mas for slice at —159° 



Note. — The 1.3 and 0.7 cm VLA-A measurements are from archival data l|Trinidad et al ] |2003l : lvan der Tak fc MentenlboOSh . The 
peak position in columns 5-6 is the relative position with respect to the reference feature of Table [2] Column 8 gives the deconvolved 
FWHM sizes of the Gaussian fits to the 1.3 cm emission and along the three axes of the 0.7 cm emission (position angles east of north). 
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Table 2. Parameters of VLBA 22.2 GHz water maser features. 



Feature 


Detection 


Vlsr 


^peak 


Ax 


Ay 


Vx 


Vy 


# 


(epochs) 


(km s-i) 


(Jy beam~i) 


(mas) 


(mas) 


(km s-i) 


(km s-l) 








North West Cluster 








1 


1 


—5.00 


14.87 


-39.72 ±0.11 


122.69 ± 0.U5 






2 


1,2,3,4 


-5.42 


14.17 


-78.51 ±0.06 


125.13 ±0.05 


4.0 ±2.9 


14.2 ±2.3 


3 


1 


—5.00 


0.55 


-46.13 ±0.08 


lz4.il lb U.UO 












Middle West Cluster 








4 


1)2,3,4 


— 11.74 


20.60 


± 0.06 


di 0.09 






5 


1,2,3,4 


—5.42 


13.41 


-4.07 ±0.12 


4. lb ± U.Ud 


4.3 ± 3.2 


14.6 ±2.1 


6 


1,2,3,4 


—6.68 


8.29 


0.02 ±0.05 


U.4/ ± U.U/ 


6.9 ± 1.9 


3.4 ± 2.2 


7 


1,2,3 


—7.10 


2.98 


-12.47 ±0.07 


U.59 ± U.Ud 


— 14.4 ± 3.3 


—0.3 ±4.1 


8 


1,2,3,4 


—4.07 


1.60 


-3.91 ±0.07 


4.I0 ± U.U7 


8.5 ± 2.1 


15.0 ± 2.2 


9 


1 


— 1.20 


0.24 


-0.37 ±0.06 


1 Tn _L c\ no 
1.79 ± U.Uo 






1 n 
iU 


A 

4 


O QO 

— z.oy 


0.24 


-16.96 ± 0.08 








11 


1 


—4.15 


0.18 


-11.78 ±0.05 


U.IU ± U.Ud 






12 


1 


—4.58 


0.17 


-3.04 ±0.05 


Q nn _L n da 
o.yy ± U.Ud 






13 


4 


—0.36 


0.14 


-14.65 ±0.07 


— U.y7 ± U.U7 






14 


4 


—6.26 


0.09 


-16.26 ±0.07 


l.Ol it U.IU 






15 


1 


— 7.94 


0.09 


-13.20 ±0.06 


1 DO _L A AQ 

1.00 ± U.Uo 






16 


1 


— 1.20 


0.06 
North 


-10.40 ± 0.07 
East Cluster 


— l.OU ± 0.06 






1 T 

17 


1,2,3,4 


—5.00 


22.25 


651.05 ± 0.06 


OArr c!r\ 1 a A£? 

305.69 ± 0.06 


—2.5 ±1.8 


4.4 ±2.1 


18 


1,2,3,4 


—9.21 


6.01 


650.18 ±0.05 


OT T AO 1 A A/^ 

317.22 ± 0.06 


—7.5 ±1.7 


16.1 ± i.U 


19 


1,2,3,4 


— 12.58 


5.44 


615.75 ±0.06 


AA C 00 1 A Ar* 

295.33 ± 0.06 


—7.5 ± 2.0 


4.3 di 2.5 


20 


1,2,3,4 


—7.53 


4.52 


618.75 ±0.10 


OAO AO 1 A 1 A 

208.02 ± 0.10 


0.3 ± 2.2 


—12.5 di 2.4 


21 


1,2,3,4 


—8.79 


2.89 


624.70 ± 0.06 


01 A OA 1 A AA 

012.39 ± U.U9 


—13.2 ± 2.1 


2.8 di 2.1 


oo 
ZZ 




—8.79 


2.29 


623.68 ± 0.11 


Q 1 Q /I _L A AK 

olz.o4 ± U.Uo 


—5.5 ± 1.7 


5.2 di 2.0 


23 


4 


—9.21 


2.03 


645.05 ± 0.06 


olD.oz ± U.Ud 






24 


1,2,3 


—8.37 


2.01 


655.31 ±0.10 


T '7 ACT 1 A A'7 

317.05 ± 0.07 


11.3 ± 4.2 


16.0 di 4.3 


25 


2,3 


—8.37 


1.85 


654.88 ± 0.08 


T T /I /I 1 A A/^ 

317.44 ± U.Ud 






OA 


1 Q /I 
1,2,3,4 


— / .00 


1.60 


608.52 ± 0.06 


Z\j ( .ID it U.UD 


— 14.0 It 1.0 


11 Q _l_ 
— il.O It Z.Z 


27=- 


1,2,3,4 


-8.37 


1.35 


620.27 ± 0.19 


311.32 ±0.16 






28 


1,2,3,4 


-7.11 


0.69 


625.15 ± 0.05 


215.20 ±0.05 


8.5 ±1.8 


-1.7 ±2.1 


29a 


1,2,3 


-6.69 


0.53 


625.44 ± 0.05 


215.62 ±0.06 






30 


1,2,3,4 


-5.00 


0.52 


648.84 ±0.08 


306.30 ±0.06 


-4.0 ± 2.0 


3.1 ±2.0 


SI'' 


1,2,3,4 


-9.21 


0.47 


641.48 ± 0.09 


315.73 ± 0.05 






32 


1,2,3,4 


-8.27 


0.42 


631.55 ±0.05 


314.48 ±0.06 


-8.5 ±1.8 


5.4 ±2.2 


33 


1,2,3,4 


-9.63 


0.40 


625.42 ± 0.09 


312.92 ±0.10 


-9.7 ±2.1 


6.4 ±2.2 


34 


1,2,3 


-8.37 


0.39 


653.85 ± 0.06 


317.29 ± 0.24 


14.1 ±3.3 


14.0 ±8.7 


35 


1,2,3,4 


-7.53 


0.39 


594.08 ± 0.05 


223.24 ± 0.13 


-14.4 ± 1.8 


8.5 ±2.2 


36 


2,3,4 


-8.37 


0.31 


657.75 ± 0.06 


309.98 ± 0.06 


24.1 ±3.0 


5.7 ±3.2 


37 


1,2,3 


-6.26 


0.24 


658.35 ±0.07 


312.14 ±0.09 


18.0 ±3.3 


-11.0 ±4.7 


38 


1,2,3,4 


-6.69 


0.23 


688.86 ± 0.08 


308.38 ±0.14 


1.4 ±1.8 


9.0 ±2.2 


39=^ 


1,2,3 


-6.69 


0.22 


685.88 ± 0.05 


311.26 ±0.05 






40 


1,2,3,4 


-6.69 


0.19 


626.85 ± 0.05 


221.95 ±0.06 


5.8 ±1.8 


-0.2 ±2.1 


41 


2,3,4 


-7.53 


0.19 


620.88 ± 0.25 


311.72 ±0.08 


-6.1 ±3.7 


0.5 ±3.0 
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Table 2 — Continued 



Feature 


Detection 


Vlsr 


Fpeak 


Ax 


Ay 


Vx 


Vy 


# 


(epochs) 


(km s ^ ) 


(Jy beam~^) 


(mas) 


(mas) 


(km s ^ ) 


(km s ^) 


42 




—7.95 


0.16 


614.48 ± 0.05 


906 1 ± 0"^ 








4 


7 


0.13 


592.54 ±0.08 


oiQ 7Q^n 1A 






44 


4 


i . oo 


0.10 


604.34 ± 0.08 


ZUU.UO ^ U.UI 












Middle East Cluster 










1 9 4 


A ^7 


3.14 


567.48 it 0.06 


cfi 97 -J- n 1 n 


y. / ^ z.z 


1 ft Q -L Q C 

— 10. y ^ 0.0 




1 9 A 


K no 

— u.uu 


0.58 


568.34 ± 0.06 


Q9 fii -1- n 1 1 
oz.uj. in u.ij. 


1 C C -L 1 Q 

j-O.u J-.y 


Q Q -L 9 7 

— 0.0 ^ Z. 1 


47 


1 9 




0.28 


565.23 ± 0.06 


1 o.uo in u.uo 


o.u ZL 0.0 










South East Cluster 










2 3 4 




4.83 


626.73 ± 0.05 


1 c:q 1 q ^ n 
lOO.-LO ZH U.UO 


7 9 4. Q 9 
— i . z zc . z 


1 . y zc o.y 


49 


9 


97 7"^ 


4.37 


621.15 ± 0.07 


1 fiA ^ n nfi 
lu^.oo m u.uu 






■^n 
ou 


]^ 


J-U. oo 


4.32 


623.54 ±0.12 


1 QR /in ^ n 07 






51 




IQ '^9 


2.05 


624.58 ± 0.08 


1 Q7 Q9 -1- 11 


7 7 A- '\ A 
— / . ( ^ 0.^ 


"lO 7 -1- S 
— ou. / ^ 0.0 


52 


2,3,4 


— 19.75 


1.71 


620.92 ±0.15 


— 1 34 25 + 1 1 


—1.9 ± 4.4 


OLI.t: — L 0.0 


53 


2 3 


— 18.90 


1.59 


621.43 ± 0.07 


j-O'i.ou u.uu 






54 


2 3 4 


9^^ 


1.52 


626.89 ± 0.06 


loz.uu ^ u.uy 


1 c; 7 _i_ Q 1 


— O'i.U ^0.0 




9 


— 17 22 


1.23 


623.17 ±0.06 


— 100. 10 ^ u.uu 








9 4. 


— 18.06 


0.97 


626.23 ± 0.08 


_-\A'i 1 -1- DQ 
— xto.-Lu ^ U.uy 






57 


9 ^ A 


—90 f^Q 
— zu.oy 


0.94 


624.09 ±0.10 


1 Q7 ^4-1-0 10 


1 n C -L Q Q 

— -Lu.u ZL o.y 


— 9 -1- 4 A 


DO 


1 9 A 


—14 27 


0.82 


612.29 ±0.23 


lie Afi -1-010 


— U.O ZL O.U 


1 n 1 -U 9 Q 
J.U. -L ^ z.y 


59 


2,3,4 


—29 44 


0.79 


624.39 ± 0.06 


1 0"^ -1- 07 
— loy.uu ^ u.u/ 


C C -|_ Q Q 
— O.U ZL 0.0 


1 A -1- ^ R 

-L'i.U ^ O.U 


60 


]^ 


—21.01 


0.70 


631.75 ± 0.14 


—AAA Oi^ n n'^ 

j.4:'±.yo ^ u.uu 








12 3 


•^n 70 


0.61 


625.12 ± 0.06 


1 c:q Oft -L n 07 
loo.yu zn u.ui 


90 7 -k R 
ZU. / O.U 


'^S A -1- A A 

00.4: ^ ^.4: 


62 


]^ 


—20.17 


0.50 


623.75 ± 0.07 


— 1 '^7 1 Q 4- n OS 
±0 ( . -Ly m U.UO 






63 


9 A 


—21.01 


0.45 


631.22 ±0.13 


— AA^ 1 9 -1- n^^ 

— XtU.-LZ ^ U.UU 


94. n _i_ K 9 

Z'i.iJ ZL U.Z 


1.6 ± 3.4 


64 


9 A 


—23.96 


0.41 


623.51 ± 0.06 


— 1 ^fi 87 -1- 07 

— J-OU.Ol ZL U.Ul 


0.4 ± 3.2 


_9c; 7 -L Q Q 
— zu. 1 ^ 0.0 


fir; 
Do 


i 


on A A 


0.38 


625.51 ± 0.06 


1 00 _i_ n c\£i 
— LOo.Jri) ± U.UD 






66 


2,3 


-15.53 


0.36 


626.12 ±0.05 


-141.52 di 0.10 






67 


1 


-22.28 


0.31 


631.30 ± 0.07 


-144.40 ±0.06 






68 


1,2,3 


-14.69 


0.24 


611.01 ± 0.13 


-114.73 ± 0.12 


-6.9 ±4.9 


13.0 ±5.4 


69 


2,3,4 


-21.85 


0.21 


630.94 ±0.05 


-144.92 ±0.05 


12.7 ±4.0 


-2.7 ±3.5 


70 


1 


-22.70 


0.20 


631.24 ±0.05 


-143.37 ±0.06 






71 


1 


-15.53 


0.18 


626.92 ± 0.06 


-143.72 ±0.07 






72 


2,3,4 


-22.27 


0.18 


630.73 ± 0.08 


-144.46 ±0.10 


11.2 ±3.8 


-3.2 ±4.3 


73 


4 


-21.01 


0.16 


623.82 ± 0.08 


-136.98 ±0.09 






74 


4 


-17.22 


0.14 


624.97 ±0.13 


-140.05 ±0.09 






75 


4 


-23.96 


0.10 


626.35 ± 0.06 


-151.82 ±0.07 






76 


4 


-25.22 


0.09 


626.43 ± 0.07 


-152.47 ±0.10 






77 


2,3 


-33.65 


0.09 


624.83 ± 0.06 


-155.15 ±0.07 






78 


1 


-13.01 


0.07 


616.11 ±0.07 


-118.00 ±0.09 






79 


4 


-22.28 


0.06 


623.72 ± 0.06 


-141.92 ± 0.10 






80 


1 


-26.91 


0.05 


625.20 ±0.20 


-159.19 ±0.65 







Reference Feature # 4: absolute position on 2009 May 6 
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Table 2 — Continued 



Feature 


Detection 


VlSR Fpeak 


Ax Ay Vx Vy 


# 


(epochs) 


(km s~^) (Jy beam"'^) 


(mas) (mas) (km s^'-) (km s""') 






R.A. (J2000) 


Decl. (J2000) 






2o'>29'"24?8228 


40°11'19'.'593 



Note. — For each identified feature belonging to a given cluster of maser emission, the label (given 
in Column 1) increases with decreasing brightness and the different epochs of detection axe reported 
in Column 2. Columns 3 and 4 report the LSR velocity and brightness of the brightest spot of eax;h 
feature, observed at the epoch underlined in Column 2. Columns 5 and 6 give the position offset at 
the first epoch of detection relative to the feature # 4 (centered on the phasc-rcfcrcncc maser channel 
# 145) in the cast and north directions, respectively. The uncertainties give the intensity-weighted 
standard deviation of the spots distribution within a feature and was combined in quadrature with 
a conservative 50 ixas uncertainty. This 50 /xas plateau accounts for deviations from a Gaussian fit 
assumption. The absolute position of the reference maser feature is reported at the bottom of the 
table and is accurate to within ±1 mas. Columns 7 and 8 report the projected components of the 
feature proper motion relative to the feature # 4, along the east and north directions, respectively. A 
proper motion of 1 mas yr~^ corresponds to 15.6 km at a distance of 3.3 kpc. 

''Features with too uncertain proper motions (> 70%) 



